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ABSTRACT 

We use numerical simulations to investigate the behaviour of the dipole moment 
of the spatial distribution of different kinds of mass tracers. We select density 
peaks of the simulated matter distribution with mean separations of 38 and 
30 Mpc to represent two samples of rich clusters of galaxies, i.e. R > 
Abell clusters and APM clusters respectively. We also extract, from the same 
simulations, samples selected to mimic the full 3D galaxy distribution of IRAS 
galaxies, and the flux-limited IRAS and QDOT galaxy samples. We compare 
the dipole moments of these "galaxy" and "cluster" samples in order to assess 
the effects of sampling uncertainties and shot-noise on the relationship between 
the "true" underlying galaxy dipole and the dipoles obtained for clusters and 
for the flux-limited galaxy samples. The results of this analysis demonstrate 
that the dipoles of both the IRAS and QDOT-like catalogues should trace 
the full 3D dipole shape fairly accurately, with the loss however of about 15- 
20% of the total 3D dipole amplitude. Furthermore, using a simple argument 
based on linear perturbation theory, on the linear biasing assumption and on 
the amplitude of the cluster dipole relative to that of galaxies, we can estimate 
their relative biasing factors quite accurately and in agreement with results 
obtained by other methods. 



Keywords: galaxies: clustering-large scale structure of Universe-dark mat- 



ter 



INTRODUCTION 



The issue we investigate in this paper is the extent to 
which different kinds of cosmic objects trace the under- 
lying mass fluctuations. We want to see how well the 
large-scale density field is traced by objects selected ac- 
cording to different criteria. In particular we will con- 
sider IRAS/QDOT galaxies and rich galaxy clusters, as 
tracers of the density field. 

It is usually assumed that number fluctuations of 
some kind of cosmic objects are related to fluctuations 
in mass M by the linear biasing factor b: 



N J 



\ M J 



(1) 



Although there are some motivations for the presence 
of a form of statistical bias relating, for example, cell- 
count variances to mean-square matter fluctuations 
through a constant b 2 (e.g. Kaiser 1984; Coles 1993), 
it should be mentioned at the outset that the motiva- 
tion for assuming a bias of the linear form (1) at each 
point is weak and this assumption may well turn out 
to be false. If velocity information is available and the 



linear biasing assumption is adopted, then it is useful to 
work with the parameter 



p = f(Qo)/b ~ nr/b, 



(2) 



where the function / is related to the rate of mass fluctu- 
ation growth in linear perturbation theory and fi is the 
present-day value of the cosmological density parame- 
ter. The method used to determine the /^-parameter, on 
which we shall concentrate here, is based on the dipole of 
the 3D distribution of extragalactic objects, its conver- 
gence and alignment properties. The 3D dipole moment 
provides an estimate of the gravitational acceleration 
acting on the Local Group (LG) and, in linear gravi- 
tational instability theory, the peculiar velocity vector, 
u(r), is proportional to, and aligned with, the peculiar 
acceleration vector, g(r): 



u(r) 



2 /(0.) g(r) _ P 
3 H Q, 4?r 



*( r )]7p dr 



(3) 



Therefore, an estimate of fi can be obtained by compar- 
ing the observed LG velocity relative to the Cosmic Mi- 
crowave Background (CMB) and the net gravitational 
force acting on it, as inferred from the observed distri- 
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bution of cosmic objects. For this method to work, the 
inferred acceleration vector must be aligned with the LG 
velocity vector. But while this is a necessary condition 
for eq.(3) to apply, it is not sufficient: one must also be 
sure that the sample of mass tracers is sufficiently deep 
that it does not miss any contribution from distant (but 
possibly very large) density fluctuations. In other words, 
the dipole vector of the distribution must converge to 
its global value within the effective depth of the sample. 

Such analyses have been performed by many 
authors, using different populations of extragalactic 
objects: optical galaxies (Lahav 1987; Plionis 1988, 
1989; Lahav, Rowan-Robinson & Lynden-Bell 1988; 
Lynden-Bell, Lahav & Burstein 1989; Hudson 1993), 
IRAS galaxies (Meiksin & Davis 1986; Yahil, Walker 
& Rowan-Robinson 1986; Villumsen & Strauss 1987; 
Strauss & Davis 1988; Yahil 1988; Rowan-Robinson 
et al. 1990; Rowan-Robinson et al. 1991; Strauss et 
al. 1992; Plionis, Coles & Catelan 1993), X-ray active 
galactic nuclei (Miyaji & Boldt 1990), X-ray clusters 
(Harmon, Lahav & Meurs 1987; Lahav et al. 1989), and 
Abell/ACO clusters (Scaramella, Vettolani & Zamorani 
1991; Plionis & Valdarnini 1991; Scaramella, Vettolani 
& Zamorani 1994; Tini Brunozzi et al. 1995; Branchini 
& Plionis 1996). The wide variety of tracers provides dif- 
ferent values of fi which can nevertheless be reconciled 
with a value of fi = 1 by invoking different levels of 
biasing with respect to the background matter distribu- 
tion. The only exception to this is the cluster dipole 
(Scaramella et al. 1991; Plionis & Valdarnini 1991), 
which, taken at face value, seems to support a value 
of Q < f. Recently, however, it has been shown that 
redshift space distortions could be responsible for this 
result (Branchini & Plionis 1996): taking such effects 
into account, a value of fi consistent with unity is pro- 
vided also by the cluster dipole (see also Plionis 1995). 

It is the purpose of this paper to investigate the 
way in which the dipole moments of different classes 
of extragalactic objects relate to that of the underlying 
galaxy distribution, and how sampling effects alter the 
convergence and alignment properties of such dipoles. In 
particular, we wish to address the following questions: 

• Do long-range contributions to the cluster dipole 
(as in the Abell/ACO sample) imply similar contribu- 
tions to the underlying galaxy dipole? 

• Does the sparse sampling of IRAS/QDOT galaxies 
seriously affect the ability of such samples to trace the 
full 3D galaxy dipole, especially when the latter has 
contributions from large scales? 

• What can be inferred from the observed dipole con- 
vergence and alignment properties of different mass- 
tracers about their relative bias parameters? 

To this end we calculate the dipole vector for vari- 
ous simulated samples constructed in order to mimic 
the flux-limited IRAS/QDOT galaxy and rich galaxy 
cluster samples. We compare these quantities with the 
dipole of the full simulated IRAS-like galaxy field in or- 
der to assess the accuracy with which the samples reflect 
the properties of the underlying "parent" distribution. 
The simulations we use are based on an optimised ver- 
sion of the Truncated Zel'dovich Approximation (TZA). 



Although this method has been demonstrated to be ex- 
tremely reliable to simulate the large-scale distribution 
of galaxy clusters (cf. Borgani et al. 1995), it cannot 
be used to predict galaxy positions with any precision. 
However, this does not represent a limitation for the fol- 
lowing analysis. In fact, the purpose of this paper is not 
to measure galaxy correlation properties, which would 
require a precise knowledge of galaxy positions, but in- 
stead to seek the effects of using different mass tracers 
to study the dynamical origin of the LG motion. It is 
sufficient, therefore, for us to use simulations such as 
these for which the large-scale matter distribution is 
accurately represented, even if the distribution of mass 
on the scale of individual galaxies is not resolved. 

The layout of this paper is as follows. In Section 2 
we set up the problem to be investigated and the meth- 
ods we will use. In Section 3 we discuss the error analy- 
sis and in Section 4 we present our results. Finally, our 
main conclusions are presented in Section 5. 



2 SETTING UP THE PROBLEM AND THE 
METHODOLOGY 

In this paper we will deal with estimates of the pecu- 
liar gravitational acceleration induced on the positions 
of suitable observers, identified in a cosmological simu- 
lation, by the distribution of matter surrounding them. 
Such estimates will be given by the dipole moment of 
the spatial distribution, around the observer, of some 
kind of mass-tracer. In order to answer the questions 
posed in the introduction we need to introduce (a) the 
moment-based method to estimate the acceleration vec- 
tor for a model simulated universe; (b) the procedure to 
identify suitable observers; (c) the method to extract 
realistic IRAS, QDOT and cluster samples for each of 
these observers. 

Note that we are not interested in studying the ef- 
fects of incomplete sky coverage, Galactic absorption, 
Virgocentric infall corrections and redshift space distor- 
tions. For such a study and its implications for dipole 
analysis, see Tini Brunozzi et al. (1995) and references 
therein. 



2.1 The formalism 

Given a set of objects at positions r;, the peculiar grav- 
itational acceleration, acting on the observer at the ori- 
gin of the coordinate system, is estimated from the 
dipole moment: 

N 

D = J2 W '*'> ( 4 ) 

8 = 1 

where i t is the unit vector pointing at the position of 
each object and N is the total number of such objects 
within the distance considered. The weights, w t , should 
be of the form w t cx r~ 2 (since both flux and gravity 
fall as r~ 2 ). 

Note that the dipole, as a function of distance from 
the observer, keeps rising till isotropy is reached, after 
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which it flattens out and converges to its final value. 
Similarly, the monopole moment M is defined as 



M 



N 

1=1 



(5) 



According to its definition, the monopole term is related 
to the number density of objects as a function of radial 
distance; in the continuous case we have that M(r) = 
47r(ra)r. Therefore, it is evident that we can define an 
estimator for the average density (n) within a distance 
R as 



M(R) 
4ttR ' 



(6) 



where R is taken to be sufficiently large (at least of the 
order of R CO nv, the dipole convergence scale), so that 
the effects of local density fluctuations cancel out in M. 
Another estimator of the mean object space density can 
be given by the inverse cube of their mean separation: 

-3 



«2 



(7) 



In the most common cases of magnitude- or flux- 
limited extragalactic object catalogues one has to take 
into account the effects of the consequent undersam- 
pling of the density field especially at large distances, 
where the radial selection functions rapidly decline. As- 
suming that the unobserved galaxies are spatially cor- 
related with those included in the catalogue, the usual 
procedure to correct for the missing population is to 
weight each observed galaxy at a distance r by l/(f>(r), 
which is the reciprocal of the portion of the luminos- 
ity function that cannot be sampled at that distance 
due to the flux or magnitude limit of the catalogue. 
Therefore the weight in eq.(4) and eq.(5) should become 
w, ~ c/> -1 (r 8 ) r~ 2 , where 



4>(n)= / <Z(L)dL (8) 

is the selection function. The lower limit of this integral, 
determined by the flux limit, is the minimum luminos- 
ity that an extragalactic object can have in order to be 
visible at a distance r, while i m ax is the maximum lumi- 
nosity of such an object. Note that the inverse selection 
function weight grows with increasing r and this acts as 
a compensating effect for the fact that poor sampling 
occurs at large redshifts. 

Finally, using linear perturbation theory (eq.3) and 
the above definition of the dipole we can relate the ob- 
served peculiar velocity of an observer, u, with that pre- 
dicted by its gravitational acceleration, v, for as long as 
the two vectors are well aligned: 



u v) = ~r~ / w > i — n = P v ( 
v ' Airn ^— ' r, 3 v 



(9) 



This relation can in principle provide an estimate of the 
/^-parameter. 

In the following analysis, we will explicitly as- 
sume that there are no appreciable density fluctua- 
tions at distances larger than our sample limiting depth 
(-Rmax = 240 ft -1 Mpc). Therefore, the density fluctu- 
ations causing the motion of the observers is assumed 



to be well sampled by the distribution of objects lying 
within i?max, which is also consistent with the observed 
convergence of the cluster dipole at ~ 170 ft -1 Mpc. 

2.2 The Simulations 

The simulations of the large-scale matter distribution 
that we use in this work are generated by using an opti- 
mized version of the Truncated Zel'dovich Approxima- 
tion (TZA) described in detail in Borgani et al. (1995). 
Here we only briefly remind the reader of its main ideas. 
Particles are moved from their initial (Lagrangian) po- 
sitions to their final (Eulerian) coordinates by the usual 
Zel'dovich mapping, x(q, t) = q + 6(t)V^>(q), where 
b(t) is the perturbation linear growth factor and ^>(q) 
is the velocity potential, which is related to the ini- 
tial density fluctuation field by the Poisson equation 
^ 72 ^( < l) = ~^(l)/ a - According to this prescription, par- 
ticles move in straight line trajectories with velocity 
v = b(t)ip(q). This gives a reasonably accurate account 
of the particle motions until "shell-crossing" occurs, i.e. 
until the mapping between initial and final positions be- 
comes degenerate. In practice, we reduce the amount of 
shell-crossing by suitably filtering the short wavelength 
modes, by using a Gaussian filter with radius chosen 
in such a way that N s = 1.1 for the average number 
of streams at each Eulerian point. Although this filter- 
ing procedure restricts the amount of non-linearity the 
TZA can cope with, nevertheless the optimized version 
of this method is well-suited to the problem of sim- 
ulating the distribution of rich clusters, in which the 
smoothing scale required is similar to the "catchment 
area" of a protocluster. 

Velocity and density fields are reconstructed on 
256 3 grid points, using as many particles, for a L = 
480 ft -1 Mpc simulation box. This size is a compro- 
mise between the need of using large samples to ensure 
dipole convergence and the need for good grid resolu- 
tion to identify mass tracers accurately. As usual, the 
box has periodic boundaries. 

The following analysis is based on a realization of a 
Cold+Hot dark matter (CHDM) model with fi = 1, 
fiftot=0.3 and fib = 0.1 for the contributions of hot 
and baryonic component respectively, H = 50 km s -1 
Mpc -1 for the Hubble constant and ag = 0.78 for the 
r.m.s. fluctuation amplitude within a top-hat sphere 
of 8/t -1 Mpc radius (see Klypin, Nolthenius & Primack 
1995, and references therein, for the cosmological rele- 
vance of this model). This spectrum normalization cor- 
responds to Qrms—ps = 20 fiK for the quadrupole of the 
CMB temperature anisotropy from the two-year COBE 
data (Gorski, Stompor & Banday 1995). 

Although different dark matter models can easily 
be simulated with our TZA method, our main inter- 
est here is not to use the dipole analysis as a test for 
such models (see Tini Brunozzi et al. 1995, for such 
an approach). This is the reason why we concentrate 
here only on the CHDM model, which we showed any- 
way to provide an excellent description of the cluster 
correlation properties (Plionis et al. 1995; Borgani et 
al. 1995). Furthermore, although many realizations of 
the same model can be quickly generated, we consider 
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here a single model realization since we verified that any 
"cosmic variance" effect on dipoles is dominated by the 
observer-to-observer scatter, rather than by differences 
between different such big simulation volumes. 



2.3 Choosing suitable observers 

We now proceed in identifying in our simulation suit- 
able, for our purpose, observers. We would ideally pre- 
fer to choose observers that are LG-like. However, as 
we shall see, we have found very few such observers and 
thus we will define observers with only some of the LG 
characteristics. The criteria for selecting LG-like "ob- 
servers" are: 

(i) Peculiar velocity of 627 ± 44 km s -1 (error cor- 
responding to 2cr uncertainties; Kogut et al. 1993) for 
a top-hat sphere of 7.5 ft -1 Mpc radius centred on the 
observer. 

(ii) Density contrast of— 0.2 < Slg < 1 within the 
same sphere. 

(iii) Dipole shape for the surrounding cluster distri- 
bution similar to the real Abell/ACO cluster dipole 
(Scaramella et al. 1991; Plionis & Valdarnini 1991; 
Branchini & Plionis 1996), i.e. a two-step increase of the 
amplitude and final dipole convergence at ~ 140 — 180 
h' 1 Mpc. 

The first two criteria are the usual ones invoked in nu- 
merical simulation studies (cf. Moscardini et al. 1995; 
Strauss et al. 1995; Tini Brunozzi et al. 1995). The last 
criterion is imposed in order to reproduce the observed 
Abell/ACO galaxy cluster dipole, since (as discussed in 
the introduction) we want to investigate (a) whether 
the particular cluster dipole shape implies that a sim- 
ilar dipole shape should be observed in the underlying 
galaxy distribution and (6) whether flux-limited IRAS 
subsamples can see such dipole contributions from large 
distances (> 120 h' 1 Mpc). 

Out of the original TV o (, s =20000 random observers 
selected in our simulation we found only 7 fulfilling all 
three criteria. Since this number is extremely small for 
a statistical study and since the questions we want to 
answer are strictly related to the last criterion, we re- 
laxed the first two, obtaining in total 106 observers. Our 
analysis will be based on this set of observers. 



2.4 Mass— tracers selection 

In what follows we present the method of extracting 
realistic distributions of clusters and galaxies from our 
simulation. Note that the distribution of simulated clus- 
ters will closely resemble the corresponding distribution 
of Abell/ACO clusters (cf. Borgani et al. 1995). How- 
ever, we will identify galaxies as simulation particles 
having the observed large-scale characteristics without 
implying that their detailed distribution and clustering 
properties should resemble those of real galaxies; our 
simulation is just not well suited for such a task. 



[tp] 

Table 1. Parameters for the Saunders et al. (1990) luminos- 
ity function. 

Parameters Values Units 

C 0.026 ± 0.008 h 3 Mpc -3 

a 1.09 ± 0.12 

(j 0.724 ± 0.031 

L. 10 (8-47±0.23) h -2 LQ 

2.4-1 Simulated Galaxy Clusters (= C) 

We identify clusters as local peaks of the evolved den- 
sity field, after reconstructing it at the grid positions. 
We select two populations of clusters; one with a mean 
intercluster separation of (dc) = 38 ft -1 Mpc (C38 here- 
after), which corresponds to the R > Abell/ACO 
sample, and one with {dc} = 30 ft -1 Mpc (C30 here- 
after), corresponding to the APM sample (Dalton et 
al. 1994). The total number of clusters within the box 
(i.e. the highest local density maxima on the grid) is 
given by Nc,t = (L I {dc)) 3 ■ Thus we have Nc,i = 2015 
and Nc,2 = 4096 for the two populations, correspond- 
ing to densities of n x 1.82 x 10 -5 h 3 Mpc -3 and 
n « 3.7 x 10 -5 h 3 Mpc -3 , respectively. 

2.4-2 Simulated IRAS Galaxies (= G) 

In order to simulate the IRAS galaxy sample we need 
to reproduce its observational flux limit. This, in turn, 
requires galaxy luminosities to be drawn from the ap- 
propriate luminosity function, 

(Saunders et al. 1990), the parameters of which are given 
in Table 1. We have used only the central values of these 
parameters in this analysis, but we do not expect the 
results to be very sensitive to changes within the known 
error limits. We should also mention at this point that 
there is possibly a significant effect due to evolution of 
the luminosity function (cf. Saunders et al. 1990). We 
shall not take this into consideration here but, if it were 
present over the redshift range we are considering here, 
it could have a big effect on the convergence properties 
of the IRAS/QDOT dipole (Plionis et al. 1993). 

To estimate the expected number density, (no), of 
the underlying distribution of IRAS-like galaxies we in- 
tegrate $(£) over the luminosity range that we want to 
represent: 

(no) = J $(i)di , (11) 

where we take i m in and i m ax to be 10 8 and 10 13 iQ 
respectively: this gives (kg) ~ 0.031 h 3 Mpc -3 . The 
lower limit in L is imposed by the fact that the very low 
luminosity galaxies are not well represented in the IRAS 
sample that we want to simulate (Rowan-Robinson et 
al. 1990, hereafter RR90). By multiplying this value by 
the volume of the simulation box used here, we obtain 
the total number of IRAS galaxies expected within the 
volume considered, No ~ 3.5 xlO 6 objects. We therefore 
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randomly select No objects from the total of 256 3 simu- 
lation particles which will represent the underlying dis- 
tribution of IRAS-like galaxies from which flux-limited 
subsamples will be extracted. Therefore, the bias factor 
of this "galaxy" population is by definition bcj = 1. 

2.4-3 Simulated flux-limited IRAS galaxies (= /) 

To simulate the IRAS flux-limited catalogue we use the 
following prescription. We take the number density of 
objects lying in a shell of width Ar at distance r to be 

N(r) = 47rr 2 (ra G )0(r)Ar, (12) 

where (ncj) is given by eq.(II) and (f>(r) is estimated 
from eq.(8) by using 

imin(r) = 47rr 2 Sii m v 60 . (13) 

Here 5i; m corresponds to the QDOT flux limit of 0.6 Jy 
and v 60 is the frequency of the passband at 60 /tin which 
is used to convert fluxes to luminosities and i m ax = 
10 13 L Q (RR90). 

The redshift distribution of the resulting Ni ~ 
12000 mock IRAS galaxies together with the theoretical 
curve of eq.(12) is shown in Figure 1. It closely resembles 
the corresponding QDOT distribution (cf. Figure 2 of 
RR90); it has a maximum at ~ 50 h' 1 Mpc (z ~ 0.018) 
and there is an extended tail out to the limiting depth 
-Rmax- The above number of IRAS galaxies expresses the 
number of objects that can be seen above the flux limit 
of 0.6 Jy. 

Evidently the number-density of objects, n, is of 
major importance in the sample definition. For this rea- 
son, we have checked the various density estimators 
eqs. (6) and (7), which give a maximum difference of 
Sn ~ 8% and are therefore in very good agreement. 

2.4-4 Simulated QDOT galaxies (= Q) 

Having generated IRAS-like catalogues down to a flux 
limit of 0.6 Jy, our next task is to build QDOT mock cat- 
alogues which constitute a random l-in-6 subsample of 
the complete catalogue. In order to check the statistical 
reliability of the results, we repeated the random l-in-6 
selection procedure 30 times for each selected observer 
using different initial seeds. Clearly, the QDOT-like cat- 
alogues have Nq ~ 2000 objects each. 

2.4.5 Assigning appropriate weights 

To estimate the cluster (C) and full IRAS (G) dipole we 
use their whole 3D distribution. Note that clusters and 
galaxies are not in fact point-like objects and thus an 
appropriate weight should really take account of their 
mass, i.e. w t ~ M 8 /r 2 . However, for sake of simplic- 
ity and uniformity, we will assume unit masses for all 
available objects. However, we verified that no serious 
discrepancy occurs in the dipole convergence, amplitude 
or alignment properties with the use of either weights. 

In the case of flux-limited IRAS samples (I&Q), we 
use the inverse selection function weights. Since it would 
be computationally expensive to solve the integral in 



eq.(8) for all the objects in each sample, we create a list 
of 240 r, distance bins having a 1 ft -1 Mpc width each 
and evaluate the integral in eq.(8) at the centre of each 
distance bin and thus recover the depth dependance of 
the selection function. Then we compare each true ob- 
ject distance rj with this list of distances and assign the 
appropriate weight to the galaxy that corresponds to n, 
as long as r t — 1 < r 3 < r t + 1. The algorithm thus allows 
us to use the previously computed and tabulated 4>(r t ) 
instead of the true <f>(rj). Clearly, this is an approxi- 
mation and it could introduce small errors. We have 
found in tests that the corresponding relative difference 
in the weights w t between the true and the above esti- 
mate of the selection function is of the order of ~ 10 -4 . 
Note that in a perfectly uniform distribution, the over- 
estimate/underestimate of c/> -1 (r 8 ) as produced by the 
above scheme, would generate similar effects, which we 
would expect to cancel out within a sufficiently large 
volume. Thus the above selection function simplification 
should not introduce any spurious modification on the 
structure, convergence and alignment of the estimated 
dipoles. 

With the above prescription to extract the different 
mass tracer populations and to define the appropriate 
weights to be used in eq.(4) and eq.(5), we proceed to 
calculate the dipole and monopole moments of the spa- 
tial distribution of these objects. To do so, we place each 
observer at the origin of the coordinate system and use 
the periodic boundaries of our simulations to ensure full 
spherical coverage. 



3 ERROR ANALYSIS 

3.1 Stability of dipole estimators 

We estimate the dipole and monopole moments for all 
the datasets extracted in the previous section using 
7.5 ft -1 Mpc, as the minimum distance, so as to be 
consistent with the definition of observer velocity and 
local density, and integrating up to 240 ft -1 Mpc (the 
simulation box limit) using bins of 10 h' 1 Mpc width. 
Note that in eq.(9) the normalization is based on the 
estimated value of the mean density of objects (rai), 
which is computed according to eq.(6) as the average 
over the last four bins, [210-240] ft -1 Mpc, to reduce 
the effects of statistical fluctuations in the value of M. 
We checked that estimating n\ over the [170-210] ft -1 
Mpc range gives very small differences in the dipole am- 
plitudes (8D ~ 1%), convergence properties and align- 
ments. 

Many authors have argued how the dipole conver- 
gence could be artificially imposed by undersampling on 
large depths which causes also the monopole to flatten. 
As can be seen from eq.(6) and eq.(9), a possible early 
leveling of the monopole moment could cause the dipole 
to flatten artificially (since the dipole gains amplitude 
much more slowly than the monopole). The monopole 
of all available datasets was found to be linearly ris- 
ing up to i?max- Note however that although this is a 
necessary condition for the reality of the dipole conver- 
gence it is not, by any means, a sufficient condition. In a 
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flux-limited catalogue you can have a linearly increasing 
monopole merely due to the volume preserving correc- 
tion factors, c/> -1 (r 8 ), although the density field may be 
extremely sparsely traced at large distances and there- 
fore the dipole could be underestimated at such dis- 
tances. As we will see there is such evidence from the 
flux-limited IRAS-like samples. 

3.2 Shot— noise estimators 

Due to the sparse sampling of the full IRAS density field 
by the flux-limited subsamples (I&Q), their moments 
are affected by discreteness effects (shot-noise errors), 
that increase with redshift because of the rapidly declin- 
ing selection function. These effects introduce a variance 
in the dipole amplitude and direction. We use two dif- 
ferent methods to estimate the shot-noise effects. 

The first method is similar to that used by Plionis 
et al. (1993). For each IRAS-like sample we generate 30 
random realizations in which we reshuffle the position 
on the sky of each galaxy while retaining its true dis- 
tance. This approach is computationally very fast and 
preserves the exact shape of IRAS selection function. 
The shot-noise dipole for each resampling is 



E 



WiTsN,i, 



(14) 



where i,sN,i are the new random directions on the sky. 

A different procedure to estimate the shot-noise 
dipole (Strauss et al. 1992; Hudson 1993; Branchini & 
Plionis 1996) is to write its components along the three 
axes, as a 



<jy, <r z . For example for the 2-axis we have 



(15) 



where a 3D = ^/ 1 <j\ + <j\ + <j 2 z and the corresponding er- 
ror along the line of sight (ID error) is a 1D = a 3D /\/3. 
By comparing the two methods we find that a 3D PS 
|Dsjv| with a maximum difference of ~ 8% for the I- 
sample and ~ 11% for the Q-sample. 

Generally we find that |Dsjv| is approximately 
(10 — 35)% of the overall I and Q sample dipole value 
while the shot-noise errors for the Q-samples are almost 
twice as large as the corresponding for the /-samples: 
this is obviously due to the much sparser sampling of 
the density field in the former case. 

Therefore, our procedure for the error analysis can 
be sketched as follows. 

• We will calculate shot-noise errors using the first 
method of eq.(14). Since it is Monte Carlo based, for 
each observer we will estimate the mean value of Dsn 
by averaging over 30 realizations of each sample. Note 
that the second method produces always shot-noise val- 
ues within ±1(7 of those provided by the first method. 

• Each estimated dipole (of the I and Q samples) will 
be corrected for the shot-noise effects by subtracting the 
ID shot-noise dipole. Therefore we will be discussing 
corrected dipole values of the form : 



• We do not apply any correction for the variance 
introduced by assuming a unit mass for all objects (cf. 
Strauss et al. 1992; Branchini & Plionis 1996). 



4 RESULTS AND DISCUSSION 

4.1 Dipole comparisons and Velocity 
fluctuations 

In order to realize a quantitative comparison between 
the different mass tracer dipole shapes, we will perform 
a correlation coefficient analysis for all possible combi- 
nations of samples. This is solely based on the shape 
of the dipole as a function of radial depth and not on 
the dipole amplitudes, which in the case of galaxies vs. 
clusters are expected to differ significantly due to their 
different bias parameters. The correlation coefficient is 
defined as 



EL<>° 



(17) 



where the quantities x t ,y t represent the different mass 
tracer dipoles and the sum is over all (or a selected range 
of) the radial bins. We will present the average value, 
(c t •), over the ensemble of the selected 106 observers. 

Furthermore, since the different galaxy dipoles 
should trace the same underlying fluctuations and have, 
by definition, bcj = hi = bq , we will assess how well the 
flux-limited samples trace the full galaxy dipole by com- 
puting their relative velocity fluctuations as a function 
of distance, averaged again over all observers: 



SVi, a (r) 



Vi(r)-V a (r) 

V G (r) 



(18) 



|D C 



|D r 



\Bsn\/V3 



(16) 



where the subscript i will indicate the I or Q samples. 
If the flux-limited samples do trace well the full galaxy 
dipole then we should find 8V(r) ~ V r. 

4-1.1 Cluster (C ) and 3D galaxy (G) dipole shapes 

We find that the dipole shapes of clusters with varying 
(dc) are extremely similar and out of the 106 LG-like 
observers only in 3 cases the corresponding dipole struc- 
tures have opposite characteristics (decreasing vs. in- 
creasing amplitudes). The correlation coefficient analy- 
sis (see Table 2) gave (c c c } ~ 0.9±0.1, while its me- 
dian value is ~ 0.94. As expected, the individual sample 
dipole amplitudes vary because of their different bias- 
ing parameters (see eq.9). Furthermore, for most of the 
106 LG-like observers the full galaxy (G) dipole closely 
matches that of the cluster dipoles. The value of the cor- 
responding correlation coefficient is (c c G ) ~ 0.7±0.3, 
with ~ 0.84 for its median. 

In Figure 2 we present a comparison between the 
G (thick line), the C38 (filled dots) and the C30 (open 
dots) dipoles for 6 characteristic observers. Note that 
the cluster dipoles are scaled downwards by a factor of 
2.5 and 2 for the C38 and C30 cases, respectively. Panels 
(a) to (d) show the most common case where the clus- 
ter dipole matches the shape of the underlying galaxy 
dipole rather well. Out of 106 observers only in 9 we do 
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detect a strong disagreement in dipole morphology, con- 
vergence and direction. Two such examples are shown 
in panels (e) and (/) of Figure 2. We have verified that 
this discrepancy should be attributed to two facts: 

• The G dipole is dominated by a very nearby, al- 
though relatively low, density enhancement which is not 
represented in the C38 clusters; 

• The galaxy dipole shape is steeply increasing within 
~ 50 h' 1 Mpc followed by a significant decrease, con- 
verging rather rapidly to a very low overall value. 

Note that in the particular case shown in panel (f ) the 
C30 dipole shape is in good agreement with that of the 
galaxy one, although the same is not true for the C38 
dipole: this just implies that the fluctuations generat- 
ing the galaxy dipole have passed the density threshold 
defining the C30 cluster sample but not that of the C38 
sample. 



4.1.2 Galaxy (G) and flux-limited (I andQ) dipole 
shapes 

Regarding the flux-limited samples (/ and Q) we obtain 
a very good correlation of their dipole shapes with that 
of the full G dipole. The correlation coefficients are al- 
ways ^ 0.85, but with a decreasing trend as a function 
of distance, which indicates that the Q and I samples 
miss some of the distant G dipole contributions. We can 
identify 3 general cases, which are (in decreasing order 
of frequency): 

(i) 6Vq,g~ t>Vi t G < 0, in ~ 66% of the observers. The 
I and Q samples underestimate the total G dipole and in 
some cases even up to ~ 30%. The dipole convergence 
depth is also underestimated, although in most cases 
the contributions to the dipole from large depths are 
apparent in both I and Q samples. 

(ii) 6Vi, a ^ SV QiG « (within 5%) in ~ 26% of the 
observers implying that in this case the I and Q samples 
trace accurately the G dipole. 

(iii) 6Vi t a — 8Vq t a > in very few cases. Most such 
observers have a very low asymptotic G-dipole value 
(;$ 200 km s -1 ) and therefore they also have a shot- 
noise dipole of the same order. 

In Figure 3 we plot the I dipole (dashed line), the Q 
dipole (dotted line) and the full G dipole (thick line) for 
the same observers as in Figure 2. In the inner panels we 
plot 8Vg,i (dashed line) and 8Vg,q (dotted line). Note 
that the observer of panel (a) is representative of case 
(ii), those of panels (b) to (d) of case (i), that of panel 
(e) of case (iii) and finally the observer of panel (/) is 
an intermediate of cases (ii) and (iii). 

In Table 2 we present an overview of the correlation 
coefficient analysis results. Deviations from the quoted 
mean values are ~ 0.25 when using all the 106 observers; 
by using the 95 observers of cases (i) and (ii) it lowers 
to ~ 0.15. Evidently, the high correlation coefficients in- 
dicate the strong correlation of dipole shapes, discussed 
previously. 

Finally we plot, in Figure 4, the 8V values as a 
function of distance for all pairs of galaxy mass tracer 



[tp] 

Table 2. Correlations coefficients, , of dipole shapes for 
pairs of different samples. Results are averaged over all 106 
observers. 





G 


I 


Q 


C*38 


C*30 


G 




0.86 


0.88 


0.69 


0.75 


I 


0.86 




0.89 


0.56 


0.63 


Q 


0.88 


0.89 




0.74 


0.76 



samples. Errorbars correspond to ±1 a estimates. A sys- 
tematic underestimate, (8V) ~ 20%, of the true dipole 
by both the I and Q samples is evident. 

4-1.3 Dipole directions and alignments 

Another important indication of how well the under- 
lying galaxy (G) dipole is traced by the dipoles of the 
different mass tracers is their respective dipole misalign- 
ment angle, A9. If linear perturbation theory of eq.(3) 
is to be used to relate velocity and acceleration, then 
a necessary prerequisite is a good alignment of the ve- 
locity and acceleration (dipole) vectors. To this end, we 
calculate the dipole misalignment angles also as a func- 
tion of distance. 

In Figure 5 (panel a) we present the A9 values as 
a function of distance for all pairs of mass tracers. As 
expected, we find that on average the distance at which 
the dipole direction converges to its final value coincides 
with of the dipole amplitude. At smaller depths, the 
misalignment angles fluctuate between (55 — 70)°. Note, 
however, that in the true cluster case the CMB and 
cluster dipoles are well aligned already from within ~ 
50 — 60 ft -1 Mpc, well before the final dipole amplitude 
convergence depth, which suggests that the LG motion 
is induced by a large-scale coherent anisotropy; see the 
discussion in Branchini & Plionis 1996 and references 
therein). 

Between the cluster samples and the G, /, Q sam- 
ples we have (A#) ~ 48°(±26°), whereas for the two 
populations of clusters (A9c 3B ,c 30 ) ~ 23°(±20°) with a 
median of ~ 19° which shows that the dipoles of the two 
cluster samples are consistent with each other in both 
shape and direction. 

The large value of A9 between clusters and galaxies 
(Fig. 5a) can be attributed in general to the sparse sam- 
pling of the underlying density field by the high peaks 
(cf. Tini Brunozzi et al. 1995) and in particular to the 
fact that the clusters, due to their large mean separa- 
tion, will tend to miss the local (near to the observer) 
contributions to the galaxy dipole. In fact if we analyse 
only those observers that trace better the local density 
fluctuations, simply requiring that their cluster dipole 
acquires its first contribution from distances ^ {d c ), we 
find (A6>) ~ 32°(±14°). 

The relevant A9 between different galaxy samples 
reveal a somewhat more complicated picture (see Fig- 
ure 5b). The most significant feature is their system- 
atic increase with radial distance which becomes more 
evident after the dipole convergence scale. We ascribe 
this to the systematic underestimate of the underly- 
ing galaxy density field due to dilute sampling, a phe- 
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nomenon which is much stronger at large distances, thus 
producing large shot-noise effects and bigger deviations. 
However, different sparse sampfes define a different de- 
gree of shot-noise amplitude, which reflects inevitably 
in the way these samples follow the 3D density field. At 
R = 180 h' 1 Mpc, we find (A6> G ,i) « 17°(±11°) and 
(A6> g ,q) « 32° (±36°), respectively. Note that the latter 
value is very close to the misalignment angle between 
the real QDOT and the CMB dipoles (RR90; Plionis 
et al. 1993). A further indication that the flux-limited 
samples underestimate the true G dipole is also provided 
by the fact that the misalignment angle increases as a 
function of distance. For example, evaluating (A#g,i) 
and (A6> g ,q) at 240 h' 1 Mpc we find ~ 26°(±19°) and 
~ 39°(±36°), respectively. 

4-1-4 General comments 

According to the above results, we can safely conclude 
that, in general, cluster samples (of Abell/ACO or APM 
type) and flux-limited galaxy samples (of IRAS/QDOT 
type) yield estimates of the dipole which are quite con- 
sistent with each other, and with the underlying dipole. 
Their amplitudes differ by a roughly constant biasing 
factor (as we will quantify below), although with a rel- 
atively large observer-to-observer scatter. However we 
have found a quite large average misalignment angle be- 
tween the cluster and galaxy dipoles. Note that in the 
observed case (Plionis & Valdarnini 1991; Scaramella 
et al. 1991; Branchini & Plionis 1996) the Abell/ACO 
cluster dipole is very well aligned with that of the CMB 
(< 25°). We find only 22 such observers in our sim- 
ulations, which constitute 21% of the whole sample, 
a rather small but non negligible fraction. For these 
observers we then find (c„ „) ~ 0.9 and (c„ .) ~ 

v c 38 > <J 1 v c 38 > 1 1 

0.8, i.e. correlation coefficients between the cluster and 
galaxy dipole shapes which are much larger than in the 
general observer case. 

In the galaxy flux-limited cases (either the I or Q 
samples) we found very good agreement of their respec- 
tive dipole shapes and alignments, although there is also 
strong evidence for an underestimation of the underly- 
ing (G) dipole contributions from large-scales. 

We summarize our conclusions on these mass trac- 
ers as follows. 

• The increase of (8V) and (A#) as a function 
of distance for the [G, /] and [G, Q] cases indicates 
that IRAS-like flux-limited samples underestimate the 
dipole contributions from large depths. 

• The value of (8Vg,i) — (SVcj,q) ~ —0.2 indicates 
that both IRAS and QDOT samples will miss roughly 
the same fractional amount of the full dipole, indepen- 
dently of the l-in-6 or 6-in-6 sampling. 

• The value of (A#) is, however, larger for the l-in-6 
than the full 6-in-6 flux-limited samples by, on average, 
~ 15°. 

4.2 A Statistical Biasing Scheme 

Various authors have used so far different methods to 
estimate the bias parameter. Correlation function anal- 



ysis (cf. Lahav, Nemiroff & Piran 1990), power spec- 
trum techniques (cf. Jing & Valdarnini 1993; Peacock & 
Dodds 1994) and ratio of count-in-cell variances at dif- 
ferent scales (cf. Borgani et al. 1995). Here we describe 
our method to estimate b which is similar to that of Plio- 
nis (1995) who estimated the cluster-IRAS relative bias 
factor using the Abell/ACO and QDOT dipoles in the 
range 10^ r^S 100 ft -1 Mpc: he found, after correcting 
for redshift space distortions, that bc,i ~ 3.5. 

As discussed in Section 2.1, assuming linear pertur- 
bation theory, we can in principle relate the simulated 
observer peculiar velocities, u, to those predicted from 
their dipole, v, as measured using eq.(9), for a particular 
class of mass tracers. Therefore, comparing observer- 
by-observer the v values for different mass tracers one 
can estimate their relative biasing. Note that since in 
the simulation we identify galaxies by a random sub- 
sampling of the dark matter (DM) particles, we have 
by definition that bcj = hi = bq = 1 for the G, /, Q bi- 
asing factors: any departure from unity value is merely 
to be ascribed to sampling fluctuations. Therefore, in 
the remaining analysis we will consider only the rela- 
tive bias between clusters and galaxies. Since fi = 1 in 
our simulation, we get from eq.(9) that 

u(r) = va(r) = b~^vc{r) , (19) 

where be is the cluster biasing factor. Therefore, we 
can obtain an estimate of be by fitting eq.(19) over the 
distance range where the dipole is estimated, by mini- 
mizing the quantity 

N , 

* ( 6 ^) = M a, ) ■ (20) 

i=i 

Here the sum is over the N radial bins used to estimate 
the cumulative mass tracer dipoles and 7 represents a 
'zero-point' uncertainty which is bound to exist due to 
the fact the mean separation is much larger for clusters 
than for galaxies and thus the cluster dipole will always 
miss the very local contributions, if such do exist (see 
Plionis 1995 for an application to real data). It is also 
evident that we cannot obtain an estimate of the good- 
ness of fit since the bins are not independent. 

Note that since eqs.(3) and (9) are of doubtful ap- 
plicability when a large dipole misalignment is mea- 
sured, we use in this analysis only those observers for 
which A#cg < 25°. We will then average the results 
over the ensemble of such observers. In Tables 3 and 
4 we list the mean values and the la scatter for the 
b and 7 parameters, respectively, as obtained from the 
^-minimization method (in which we use a t = 1). The 
average formal fitting error in b and 7 was found to be 
~ 0.04 and ~ 200 km s -1 , respectively. 

In order to check the reliability of using the dipole 
amplitudes to estimate the cluster biasing factor, we 
computed it also by the count-in-cell method (see also 
Tini Brunozzi et al. 1995). We randomly displaced 5000 
spheres within the simulation box, with radii rang- 
ing from 50/t -1 Mpc to 120/t -1 Mpc and estimated the 
r.m.s. value, a c (R) of the cluster count within such 
spheres. Therefore, the cluster biasing factor can be sta- 
tistically defined as 
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Table 3. Relative biases between clusters and galaxies. 



SAMPLES 

C*38 
C*30 



G I Q 

4.1± 1.1 4.2± 1.6 4.8± 1.6 

3.3± 0.6 3.4± 0.9 3.8± 0.8 

[tp] 

Table 4. Relative 7 parameters (in km s — 1 ) between clusters 
and galaxies. 

SAMPLES G I Q 

C 3S -520±360 -290±500 -460±315 

C 30 -220±390 -50±440 -170±380 



(21) 



Here cr DM is the r.m.s. fluctuation for the DM distri- 
bution, which we computed from the power spectrum 
P(k) according to 



-!- / dkk 2 P(k)W 2 (kR) 
2tt 2 ' 



1/2 



(22) 



with W(kR) the Fourier transform of the top-hat 
sphere, W(x) = 3(sin kR - kR cos kR)/(kR) 3 . We find 
that the cluster biasing factor is constant to a good 
accuracy over the whole scale range, having a value 
bc 3fs = 3.9 ± 0.2, in nice agreement with the result re- 
ported in Table 3. 

A further point in this analysis concerns the large 
values of 7 (Table 4) which one would have naively 
expected to be rather close to zero. However as we 
discussed earlier, these 7 values just reflect the clus- 
ter dipole 'zero-point' uncertainty which is due to the 
fact the cluster distribution cannot in general accurately 
trace the galaxy distribution near the observer, there- 
fore missing local dipole contributions. A similar situa- 
tion appears in the real Abell/ACO cluster case, where 
the cluster dipole cannot trace the contribution to the 
LG motion of the Virgo cluster; see Plionis 1995 for an 
attempt to correct for such effects. 

However, the method we are using is only a sta- 
tistical measure of the relative biases, based on vari- 
ous assumptions posed by the corresponding numerical 
experiments. Although the various uncertainties (lumi- 
nosity function variations, large shot-noise errors due 
to sparse sampling, etc.) may cause a possible poor es- 
timate of the relative biases, when added together, the 
fact that we actually manage to recover the mean bias 
factors so close to those measured by other methods, 
should be considered a success for this approach. 



5 CONCLUSIONS 

We have used numerical simulations of a Cold+Hot dark 
matter model, based on the Truncated Zel'dovich Ap- 
proximation, to simulate the distribution of rich clus- 
ters, of IRAS galaxies and of flux-limited subsamples of 
IRAS galaxies in order to examine how well the different 
mass tracers reflect the underlying "galaxy" dipole. We 
address also the effects of the QDOT sparse sampling 



strategy to the determination of the underlying dipole 
structure. 

Two cluster samples were generated; the first with 
a mean intercluster separation of {dc} = 38 ft -1 Mpc 
(in order to mimic Abell/ACO sample) and the second 
with (d c ) = 30 h' 1 Mpc (resembling the APM sam- 
ple). Galaxies are identified from a random subsampling 
of the whole DM particle distribution, having, however, 
the same average density and radial selection functions 
as the real IRAS and QDOT samples. We have allowed 
ourselves such a simplification of the galaxy identifica- 
tion procedure since we are interested only in analysing 
differences between the dipole structures of different 
mass-tracer distributions, rather than in the details of 
the distribution of "true" simulated galaxies in a spe- 
cific DM model. We chose observers measuring a cluster 
dipole similar to the Abell/ACO one; then we calcu- 
lated for each of them the dipoles for all the catalogues 
of mass tracers. Care has also been taken to estimate 
and correct for the relative shot-noise errors generated 
by the sparse sampling. 

Using a a statistical correlation coefficient analy- 
sis we found that in most cases the cluster dipole shape 
reflects a similar 3D galaxy dipole shape while the align- 
ment of the two vectors is fair but not extremely good, 
(A#c 38 ,g) ~ 48°. Such large misalignment angles should 
be attributed to the fact that clusters are high peaks of 
the underlying density field and therefore, due to their 
large mean separation, they inevitably tend to miss the 
local (i.e. near to the observer) gravitational contribu- 
tions, which in most cases dominate the galaxy dipole. 
Consistently, if we restrict the analysis to those ob- 
servers for which the cluster distribution traces fairly 
well the local galaxy dipole contributions we systemat- 
ically obtain a better alignment, (A#c 38 ,g) ~ 32°. 

Comparing the dipoles obtained from the flux- 
limited galaxy samples and from the whole galaxy dis- 
tribution, we find a close match of their shapes for 
both IRAS- and QDOT-like samples. However, both 
the flux-limited samples miss a fraction of the total 
dipole amplitude (on average ~ 15% — 20%). The aver- 
age misalignment angle between the dipoles of the whole 
distribution and of the IRAS-like sample is (A#) ~ 17°; 
for the QDOT-like case it is (A#) ~ 32°, similar to 
the observed value between QDOT and CMB dipoles 
(RR90; Plionis et al. 1993). It should be noted that the 
quoted misalignment angles between the galaxy samples 
are not shot-noise-free; such corrections have only been 
applied to their dipole amplitude estimates. 

As a main conclusion of this paper, we predict that 
the dipole of the full 6-in-6 QDOT sample will have 
about the same amplitude as the original l-in-6 sample, 
or in other words it will underestimate the true dipole 
by the same amount as the l-in-6 sample (by ~ 15% — 
20%). However, the alignment of the 6-in-6 dipole with 
that of the CMB should become better and within ~ 
10° - 17°. 

Finally, we have used a statistical method, based 
on linear perturbation theory and on the linear biasing 
assumption given by eq.(l), to estimate the relative bias 
factors of the different mass tracers. Our results agree 
well with those based on the count-in-cell statistics, in- 
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dicating that our method is robust, although it produces 
quite a large observer-to-observer scatter in the derived 
biasing parameters. 
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FIGURE CAPTIONS 

Figure 1. Number density of the mock IRAS galaxies 
as a function of distance. The superimposed line cor- 
responds to that expected in a perfectly homogeneous 
universe. 

Figure 2. Cumulative dipole amplitudes for 6 different 
observers. The C38 and C30 cluster dipoles are repre- 
sented by the filled and open dots, respectively. The 
thick line corresponds to the full underline galaxy (G) 
dipole. 

Figure 3. The G (thick line), I (dotted line) and Q 
(dashed line) galaxy dipoles for the same observers as 
in Figure 2. The inner panels show their velocity fluc- 
tuations SV. 

Figure 4. Mean velocity fluctuations, (SV), as a func- 
tion of radial depth. Errorbars are la scatter (shown 
only for a few points for clarity). 

Figure 5. Mean misalignment angles, (A#), between 
C38 clusters and galaxies (panel a), and between the 
galactic samples, G vs. I and G vs. I (panel b), are 
shown as a function of depth. The squares in panel (a) 
represent (A#) between the C38 and C30 clusters. Er- 
rorbars, plotted only for one set of points, are similar in 
the other cases. 



